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Abstract: We compare the amplitudes for the long-distance scattering of three gravi- 
tons in eleven dimensional supergravity and matrix theory at finite A^. We show that 
the leading supergravity term arises from loop contributions to the matrix theory ef- 
fective action that are not required to vanish by supersymmetry. We evaluate in detail 
one type of diagram — the setting sun with only massive propagators — reproducing the 
supergravity behavior. 
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1. Introduction 

In the original work of BFSS [|l|, it was conjectured that M theory could be described 
by a matrix model in the large limit; later, refs. [0, ^, Q gave meaning to the matrix 
model also at finite as describing M theory on the discrete light-cone. Even though, 
in the long-distance regime, M theory is supergravity in eleven dimensions, there is no 
guarantee that the supergravity result will always match matrix theory at finite N due 
to the presence of a small scale in the problem [^. For this reason, it is of importance 
to test in actual computations how far the two agree. 

The most convincing test to date has been the comparison of the two-body scat- 
tering @, 1^ . ^ An even more stringent test would come from multi-particle scattering. 
The case of three 11-dimensional gravitons, carrying Kaluza-Klein momentum in the 



10 compactified direction, has been considered in [11 1 where the authors claim that 



a term present in the supergravity amplitude cannot arise in the matrix model. 

In this note we discuss this important issue by reconsidering the computation of [|ll 
Our result is encouraging for matrix theory: contrarily to what reported, we find that 
there are (O-l-l)-dimensional Yang-Mills (YM) graphs which lead to the same behavior 
with respect to the relative distances and relative velocities as in the supergravity result. 

The diagrams we have considered are two-loop diagrams in the bosonic sector — there 
are various similar diagrams which can give rise to the same behavior — and we have 
analyzed in detail one of them, the setting-sun diagram with all massive propagators, 
which only arises in the three-body problem. 

When considered in the framework of the effective action arising at one loop by 
integrating out the long-distance degrees of freedom (the heavy modes), our result 



^For earlier works on different aspects of this problem, see |T^. 
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originates from a term that does not vanish after summing over all bosonic and fermionic 
contributions — even though the effective operator for the remaining light modes does 
not contain any explicit dependence on the velocity. In the final two-loop effective 
action, obtained after integrating out all modes, the velocity independent terms cancel 
(in agreement with supersymmetry |^2|, |13|) while, the relevant contribution arises from 
a term of order which is not expected to vanish by supersymmetry. 

We have not attempted to compute the numerical coefficient, which would require 
the algebraic sum of the various bosonic, fermionic and ghost diagrams of the YM theory 
at order v^. Therefore we present here a minimal result, which we feel nonetheless to 
be important because of the recent discussion concerning whether one could or could 
not find the supergravity behavior in the diagrams of the YM formulation of matrix 
theory at finite N. 



2. The amplitude in supergravity 

The simplest way to obtain supergravity amplitudes is by means of string theory. Since 
it is a tree-level amplitude, it is consistent with conformal invariance in any dimension- 
ality, in particular in D = 11. We consider the bona fide superstring theory (where there 
is no tachyon) and the scattering amplitude of three (11-dimensional) gravitons, and 
look at suitable pinching limits, where only intermediate massless states are coupled 
to the external gravitons. Those states are themselves 11-dimensional gravitons. We 
then compactify the 10**^ space dimension giving mass to the external gravitons, which 
will thus correspond to 10-dimensional DO-branes. Keeping zero momentum transfer 
in the 10**^ dimension, the intermediate states remain massless and correspond to the 
various massless fields of 10-dimensional supergravity. 

The supergravity amplitude is thus obtained from that of superstring theory by a 
limiting procedure that isolates the relevant corners of the moduli space. We follow 
1^], where the appropriate technology is explicitly developed. 

By considering only the part of the complete amplitude that is proportional to 



£1 -e'lea -4^3 (2-1) 

e being the external graviton polarization tensor, we obtain the amplitude Aq for six 
graviton vertices: 
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The eleven-dimensional momenta are chosen to be 



= (E,, p, - q,/2, M,) p'^ = {~E[, -p, - q,/2, -M,) (2.3) 

where pf = 0, Ei ~ Mj + (pj — qj/2)^/2Mj and Mj = Ni/Ru are the momenta 
in the compactified dimension. Energy-momentum conservation gives I]i = and 

Ei Pi ■ qi = 0. 

In order to obtain a non- vanishing result in the field theory limit [a' — > 0), we must 
extract three poles in momenta, each of them bringing down one power of («')~^ and 
thus compensating for the three powers of a' in front of A^; each pole originates from 
a pinching limit in which some of the Koba-Nielsen variables come close to each other. 
The pinching limits corresponding to the grouping of the six external vertices into three 
pairs with i = 1, 2 and 3, give field theory diagrams where the incoming (i) and 
outgoing [i') particles of each pair describe the world-line of the DO-brane number i. 

In particular, we are interested in seven pinching limits. One of them corresponds 
to the so-called Y diagram where each of the three world-lines are coupled to one 
intermediate massless state, and the three intermediate states meet at a point. We 
call the corresponding amplitude Ay- In addition, there are the diagrams where one 
world-line is coupled to two intermediate states, each of them attached to one of the two 
other world-lines (with six possible choices). We disregard terms, which are interpreted 
as re-scattering effects, where a world-line interacts successively with two intermediate 
states, with an external particle's propagator in between. Thus, besides the Y diagram, 
we are left with diagrams where two intermediate states originate from the same point 
of one of the world-lines. We denote the corresponding amplitude by Ay and, therefore, 
we have that Aq = Ay + We keep only those terms giving the maximal singularity 
in the momentum transfers. 

Let us first consider the amplitude (taking, for the moment, A^j = 1) 



A^J = 2 El ■ e[ 62 ■ 4 ^3 ■ ^3 
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Eq. (|2.4|) is the same supergravity amplitude considered in [|TI 



Eq. (|2.4]) contains three possible singular configurations in which two of the q^'s are 
small, describing by Fourier transform the large- distance interaction of one DO-brane 
with the other two. Consider the case where the distance of the brane number 1 from 
the branes number 2 or 3 is much larger than the distance between the branes 2 and 
3. The two singular terms proportional to 

or (2.5) 

qiq2 qfqs 

are relevant. Let us consider the first one, as the second one is obtained by interchanging 
2 with 3. The Fourier transform gives (rj being the position in space of the i-th brane) 

ay = 2 El- e[ 82 ■ 4 ^3 ■ 4 (P3 - P2)' (Pa - Pi)' (P2 - Pi)' 
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To get the case of generic Ni, we have to replace 



(P.-P,)^-^P? + ^P?-2P..P, (2,7) 

and write the momenta in terms of the velocities as pj = MjVj while bearing in mind 
that Mi Ni. We normalize the amplitude by dividing the result by the product of 
the Mi and obtain: 

N^N2N3{^r^ - V2)2(v3 - Vi)2(v2 - Vi)^ 

civ ~ , TTi • (2.8j 

|ri - raTlra - r^]' 

In order to compare ( p.8| ) with matrix theory we consider the eikonal expression 
where we integrate over the time t along the world-line trajectories. For simplicity 
we take the velocities of all three particles to be along the axis and the relative 
displacements to be purely transverse (impact parameters). In other words, for the 
i-th particle, = [viUit + b,;), with bj • fii = 0. This integral gives, in the limit where 
5 = |bi - bal > 6 = |b2 - bg], 

N,N2Nsvl,vlv!, N^N2Ns\v23\vhv!2 

where Vij = Vi — vj. The other term in ( p.5| ) gives the same amplitude with B replaced 
hy B' = |bi — bsl ~ 5. It is the amplitude civ in ( p. 91) that we want to reproduce in 
the matrix theory computation. 

As for the Y diagram, the corresponding eikonal expression ay turns out to be 
sub-leading in our limit (see the appendix A). 

3. The amplitude in matrix theory 

The derivation of the Feynman rules and the computation of the relevant diagrams 
follow closely those of M. We use units where 



5'YM 



{Ru/\lf" = l, (3.1) 



the quantities Ru, Ap and qyu being the compactification radius, the Planck length 
and the Yang-Mills coupling, respectively. The bosonic part of the gauge fixed action 
reads 

S = J dt Tr ^(ig + i; • + 4 i i?fc [ao, Xk] - [Rk, a-of - [Rk, Xjf 
+2 i Xk [oo, Xk] + 2 oo] [ao, Xk] - 2 Xj] [x^, Xj] 
-\aQ,Xkf - -[xk.Xjf^, (3.2) 
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where and Xk are hermitian matrices representing the fluctuations and Rk is the 
background. Since we are studying the scattering of three DO-branes, with two in- 
dependent velocities and impact parameters, we need to consider at least a rank two 
group, namely SU{3). The overall factors of Ni, representing the longitudinal momen- 
tum will be fixed at the end. 

We choose the same background that led to (p^), namely, 

/vit \ \ 

Ri = I V2t \ and Rk = \ bl \ k > 1. (3.3) 
VO v^tj \0 bl) 

We factor out the motion of the center of mass by imposing vi + V2 + vs = and 
bl + b2 + b3 = 0. 

We use a Cartan basis for SU{3), where and denote the generators of the 
Cartan sub-algebra and (a = ia^, ±a^, ±a^) the roots. We also define the space 
vectors 

R° = ^ a"Tr (i/'^R). (3.4) 

a=l,2 

With the standard choice of H"^ and a, this definition singles out the relative velocities 
and impact parameters, e.g. i?"^ = (^'2 — ^3)^ = f^^t and, for k > 1, R'f = b\ — b\ = b"^ 
together with cyclic permutations. According to the previous section we choose the 
relative distance of the first particle with the other two to be much larger than the 
relative distance of particle two and three, in other words, we set 

|b°'| ^ |b"'| ^ B > |b"'| ^ 6 and Bb:^v, (3.5) 

where in our units v has the same dimensions as 6^. 

The propagators and vertices can be easily worked out from the gauge fixed action 
( |3.2| ), with two points worth stressing: first, the quadratic part (yielding the propa- 
gators) is diagonal in root space; second, contrary to the SU (2) case, there are now 
vertices with three massive particles (corresponding to the three different roots). The 
second point is particularly crucial because it is from a diagram containing those ver- 
tices that we find the supergravity term. 

We find twenty real massless bosons and thirty massive complex bosons. We only 
need consider some of the latter to construct the diagram. Writing Xk = x^H"' + x^Ea, 
with = x^*, we define the propagators as 

{xrih)xnt2)) = A''{ti,t2)kr (3.6) 

The fluctuation xi associated to the background Ri mixes with the field oq (the fluc- 
tuation of the gauge potential). 

We focus on the vertex contained in the term of the effective action (|3.2|) of type 



-2 Tr ([Ri,Xj][xi,Xjfj , 



(3.7) 



which gives a vertex with two massive bosons and a massless one and another one 
with all three massive bosons. Focusing on the second case and choosing a particular 
combination of the roots we obtain a term of the type 

■ ■ x"' = v^' t xl^xfxf + hfx]\fxf , (3.8) 

The two-loop setting-sun diagram in Fig. 1, which is obtained from the insertion of 

b 
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Figure 1: The settmg-sun diagram, 
two of the vertices ( |3.8| ), yields a contribution given, up to overall numerical factors, 

by 

ae = jdT jdt (R"^(ti) ■ x"'*(ti) ■ x-'(t2))(x"^*(ti) ■ x-'*(ti)x"'(t2) ■ x"'(t2)) 

(3.9) 

In eq. (gj), we defined t = {ti - t^^jl, T = {t^ + 12)/2. 

Our strategy is to do first the integration over the time difference t. Because dis- 
tances are large, this leaves us with a local expression C in the R"(T) and their deriva- 
tives. 

The propagator ( |3.(j| ) can be written as 

(A")h = — / ^= WUu; e-^*' , (3.10) 

rt"^ Jo y/A-KU 

where Wki is defined in the appendix B. 

We note first that the integration over u can be extended to infinity, up to expo- 
nentially negligible corrections for (i?")^/i?" oo. 

2 3 

The important point is that for the heavy propagators A" and A" , we can put 
Wki = Ski, up to terms of order not less than {R°'^''^Y / [R^^'^Y which are not relevant 
for our computation, as it will become clear in the following. Thus, to the order in 
which we are interested, we can replace the heavy propagators with those obtained as 

2 3 

if i?" and were constant, that is 

(A"'V = '5./^e-2^"'l*l, (3.11) 



and similarly for A"^. Further, with the same accuracy, we can let ~ {12) ~ 

R^^{T), since t < + R"^), and, for the same reason we can put A'^^{T,t) ~ 

A"'(T,t = 0). 

We can now perform the integration over t, obtaining 

C ^ ^'3 2 ^ 3(x"'*(T) -x^Vt)). (3.12) 

where i?" still depend on T. 

A comment is in order. The effective lagrangian (|3.12| ) looks like a one-loop correc- 
tion to the velocity independent local effective potential: 

5V,-ioop = F [i?"(T)] x"'*(T) ■ x"'(T) . (3.13) 

One could ask whether the potential (|3.13|) should vanish by supersymmetry rea- 
sons, once one performs the sum over all the one-loop diagrams made of heavy fields 
(including the gauge, ghost and fermion fields). The answer is: no, it is not canceled. 
In order to verify that this is not the case, we have evaluated the velocity independent 
part of the sum of all the diagrams with a;"^*x"^ external lines, for a generic constant 
background, and found that the term of eq. ( p.l2| ), coming from the heavy loop of our 
setting-sun diagram, does not cancel. In fact, the other heavy loop diagrams carry a 
different dependence on i?". This computation is easily done in configuration space, 
because of the simple expression of the velocity independent propagators, see eq. ( |3.11| ). 

As a further check, we have also evaluated the velocity-independent effective po- 
tential at two loops as a function of the background — that is after performing the 
integration over the heavy and light fields as well — and found it to be zero, as expected 
by supersymmetry. Evidently, the cancellations due to the symmetry of the theory do 
not occur at the level of the first step of integrating only over the heavy fields. 

Coming back to our computation, the last propag atorA" in (|3l2|) can be evaluated 
by expanding W in powers of u in eq. ( |3.1CI|) . The relevant terms are those of the kind 
(recall the overall factor 1 / R"^ in front of (|3.10|) ) : 

Performing the last integral over u and substituting into ( p.l2|) we obtain 

A = JdTC = JdT —^—^^—^——^^ i—^ . (3.15) 

Notice that A is proportional to the sixth power of R and therefore is the first term in 
the expansion in velocities for which no obvious supersymmetry-induced cancellation 
is expected |]r2|, |T3|. 



To verify that the action ( p.l5|) contains (|2.9| ) one selects terms proportional to 

[R"'^ R^^Y = "^13 "^12! ^"^^ performs the integration over T to get a term of the form 

A^'"J^lATr , „,, ~%|^. (3,16) 



The appropriate powers of A^, can be deduced — following [0 — from the double-line 
notation in which the setting-sun diagram is of order A^^; this factor must be N1N2N2 
for the diagram to involve all three particles. We thus find the term 

® w ^ ' 



which reproduces the behavior of the supergravity result ( p.9|) , that is, clq ~ Ov- 

Of course, to verify that matrix theory really matches supergravity, one should also 
check the numerical coefficient as well as the matching of other terms with different 
powers of B, b and Vij. That would require considering all the various graphs of the 
YM theory. 
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A. The Y diagram 

The Y diagram gives a term in the amplitude of the from 

A ^ A I I / I 



X UP2 - PsJ 



+ (Pl - P2) 



' q? qi qi 






- Pl) + q2 • 


(pi- 


- P2) 


- P3) + qi 


(pi 


-P2) 


- Ps) + qi 


(P3 


-Pl) 



(A.l) 



Notice that Ay = whenever = pj, as it is also true for As/. 

In the limiting regime (|3.5|) , and in our particular kinematic configuration, the term 
( |A.1|) above turns out to give a contribution proportional to 



ay 



J dt ndq, e^^- (^--*+^^)5(Eq,)Ay = /(v) , (A.2) 



where f{v) is a homogeneous function of degree five in the velocities Vij, and it is 
therefore sub-leading with respect to in (^.91). 



B. The propagators 

The explicit form of the propagator is 



^kl 



27rsinh2f"s 
2 sinh^ f 



exp I — f " (t^ coth v"s + tanh f "s) | 
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(B.l) 



By changing the integration variable to 

sinhw^s 



u = 



i?" coshu°s 

we can rewrite the propagator (|B.1|) as in eq. (|3.10|) , with 



(B.2) 
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The important point about W is that it can be expanded in powers of 



(B.3) 
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